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Abstract: DIS data from HERA show a striking regularity as a 1 p is a function of the ratio 
r = Q 2 /Q 2 (x) only. The scaling function shows a break at r ~ 1, which has been taken 
as an indication for saturation. However, besides saturation also the transition between 
dominance of fej_-ordered (DGLAP) and £;j_-non-ordered (BFKL) evolution contributes to 
a break around this value of r, as well as the suppression for small Q 2 due to finite quark 
masses and confinement. In this paper we use a dipole cascade model based on Mueller's 
dipole model, which also includes energy conservation and pomeron mergings, to investigate 
the contributions of these different effects to the scaling behaviour. As a result we predict 
that the scaling function for r < 1 will be modified when data for Q 2 > 1 GeV 2 become 
available. We also investigate the scaling properties of the charm contribution and the 
impact parameter dependence of the saturation scale. 
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1. Introduction 

Data from deep inelastic scattering (DIS) experiments at small- 2; exhibits an interesting 
property called geometric scaling [1]. This means that the total 7*p cross section is not a 
function of the two variables x and Q 2 separately but rather a function of the combination 
Q 2 /Q 2 {x) only, where the "saturation scale" Q s is defined such that saturation is expected 
to occur at Q- values below Q s . 

DIS is quite successfully described by the Golec-Biernat-Wiisthoff (GBW) model [2,3] 
in which the virtual photon splits into a qq dipole long before the interaction with the 
proton. The GBW model is also called the saturation model, since it explicitely assumes 
that the dipole-proton cross section, a^p, saturates to a constant value do as the dipole 
size, r, gets large. To be more precise, the GBW model assumes that 

<Td P = <t {1 - exp(-r 2 /AR 2 (x))}, (1.1) 
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where the "saturation radius" Rq(x), identified with Q s l {x), decreases with decreasing x. 
In the GBW fit Q 2 has the form 

Ro 2 = Qlcw = Qo with Q = lGeV, x = 3-10- 4 , A = 0.29. (1.2) 

We thus see that cidp is a function of r/Ro(x) only, and consequently it satisfies "geometric 
scaling" . For massless quarks this also implies that the 7*p cross section is a scaling function 
of r = Q 2 R 2 = Q 2 /Ql GW , and this feature is indeed confirmed by experimental data as 
demonstrated in ref. [1]. (This reference also finds scaling if the relation lni? 2 , ~ Alnx is 
replaced by lni? 2 , ~ — x~ 08 .) 

Plotting the observed j*p cross section as a function of t one can also see that the 
growth of the cross section towards smaller r is reduced noticeably when r < 1, correspond- 
ing to values satisfying Q 2 < Q 2 GW . Later analyses including DGLAP evolution for larger 
Q 2 [4] improve the agreement with the HERA data, with the consequence that saturation 
is beginning at somewhat smaller r-values. This good agreement together with the success 
of the GBW model in describing diffractive data, has been taken by some authors as a 
proof that saturation exists, and that it has been observed at HERA. 

Including the finite quark mass in the wavefunction describing the 7* — qq coupling 
introduces scale-breaking effects. In particular the charm quark gives a large contribution 
to the cross section, which is phased out for Q 2 < Am 2 . In ref. [2] a fit including the charm 
quark mass actually shifts the expected onset of saturation to much smaller x-values, with 
xo = 0.4- 10 -4 in the expression for Rq or Q s in eq. ( |1.2|) . Such a shift is confirmed in 
the analysis in ref. [5], which includes the mass of the charm and beauty quarks into the 
DGLAP-improved formalism of ref. [4]. 

For small Q 2 also the mass of the light quarks becomes important. In refs. [1,2] this is 
taken into account by replacing the variable x = Q 2 /W 2 in the definition of Ro or Q s ,gw 
by 

. , Q 2 + 4m 2 f 
x = x(l + 4m 2 /Q 2 )= w2 f . (1.3) 

In this way the scaling relation can be studied also for very small Q 2 and for photoproduc- 
tion. 

The test of scaling in the region of small r is, however, limited by the fact that small 
values of r are reached only for small Q 2 , since larger Q 2 -values would need energies not 
accessible at the HERA accelerator. Thus the data for r < 0.5 — 1, where the change 
in the slope of the j*p cross section is observed, are all obtained for Q 2 -values smaller 
than 1 GeV 2 , which means that they are all in the non-perturbative region. The limited 
kinematical range at HERA also implies that there is little overlap between data at different 
Q 2 for fixed r, which implies that it is relatively easy to achieve a scaling result by adjusting 
the quark masses. 

The question of scaling for Q 2 > Q 2 is discussed by Iancu et al. in ref. [6]. These 
authors argue that if Q 2 is defined as the scale where the scattering probability is of order 
1, then the BFKL evolution equation implies that the quark and gluon distributions have 
to obey geometric scaling in the range 1 < \n(Q 2 /Q 2 ) -C \n(Q 2 / f& CD ). This rather wide 
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range in Q 2 results from the fast diffusion in In kj_ and the fast growth towards small x in 
the leading log BFKL evolution. Beyond this range in Q 2 the BFKL diffusion is gradually 
replaced by the (not explicitely scaling) double leading log result. 

A numerical analysis of the diffusion in the BK equation is presented by Golec-Biernat 
et al. [7], including non-leading effects from a running coupling and from the so called 
kinematical constraint. The result from a running a s is that lnQ^ for high energies grows 
~ y/Y, rather than proportional to Y, as assumed in eq. ( |1.2| ). It also significantly reduces 
the diffusion into the region of large k\ . The kinematical constraint reduces this diffusion 
further, and these results therefore put a question mark for the large scaling region domi- 
nated by BFKL dynamics obtained in ref. [6]. A related work is presented by Kwiecihski 
and Stasto [8], where they study DGLAP evolution starting from scaling initial conditions 
on the line Q 2 = Q 2 , with Q 2 determined by some unspecified dynamics. The result from 
this approach is also that scaling is approximately preserved in a large domain above the 
line Q 2 = Q 2 S . 

In the past few years it has been observed that a scaling feature is inherent in the 
asymptotic solutions to the evolution equations in high energy QCD. It was also realized 
that the non-linear Balitsky-Kovchegov (BK) equation, which is the mean field version of 
the more general Balitsky-JIMWLK hierarchy (B-JIMWLK), is similar to a certain type 
of equation, well known in statistical physics, called the Fisher-Kolmogorov-Petrovsky- 
Piscounov (FKPP) equation, which is known to have traveling wave solutions [9,10]. Writ- 
ten for a function u(x,t), which depends on x and the time t, the solution for large t has 
the form of a traveling wave, u(x — vt), where v is the speed of the wave. The similarity is 
expressed by the fact that the BK equation lies in the same universality class as the FKPP 
equation. 

More recently, the importance of fluctuations in small- x evolution has been better 
understood [11], and this has led to the modification of the original version of the B- 
JIMWLK hierarchy into a new hierarchy of equations, which include both fluctuations and 
saturation effects. These new equations are also referred to as pomeron loop equations, 
since they contain both pomeron splittings and pomeron mergings in the evolution. This 
hierarchy of equations can also be written as a single Langevin equation, which is very 
similar to what is called the stochastic-FKPP (s-FKPP) equation [11]. The study of the 
asymptotic behaviour of the solution to this equation leads to the prediction of a new type 
of scaling law, called diffusive scaling [12], which is expected to hold at very high energies. 
It ought to be emphasized that the traveling wave solutions discussed here are asymptotic 
solutions expected to be relevant at extremely high energies, and therefore cannot be used 
to explain the scaling observed in the HERA energy regime. 

The above discussion raises three important questions: 

i) What is the importance of saturation for the scaling behaviour? 

ii) What is the dynamical mechanism behind the scaling observed for r > 1, i.e. for 

Q 2 > 

iii) Will the cross section still be scaling for r < 1, when data for larger Q 2 are 
available? 
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In this paper we will argue that geometric scaling is expected also in the absence of 
saturation, and not only in the region dominated by BFKL diffusion, but also in the dou- 
ble leading log domain, where A;j_-ordered (DGLAP) evolution chains are most important. 
Scaling appears naturally in a dipole cascade model [13,14], which is based on Mueller's 
dipole evolution [15-17] but also includes energy-momentum conservation, pomeron merg- 
ing, and a simple model for the proton. The MC implementation of the model reproduces 
both F2 data from HERA and the total cross section in proton-proton scattering. This 
model shows geometric scaling for Q 2 below as well as above Q s GW , and for the one 
pomeron contribution as well as for the full unitarized result. In the model the transition 
between BFKL diffusion and fcj_-ordered evolution occurs for Q 2 quite close to Q 2 Gw , and 
there are three different effects which all contribute to the change in the scaling curve for 
r ~ 0.5 — 1: Saturation, the BFKL-DGLAP transition, and the finite effective masses for 
the light u-, d-, and s-quarks. 

In order to get an intuitive understanding of the scaling feature we will besides the 
results of the MC simulations also discuss two simple approximations, which contain the 
basic features of DGLAP-BFKL evolution and the colour dipole cascades. 

Diffractive excitation or rapidity gap events correspond to a large fraction of the events 
at HERA. Diffractive scattering is related to the impact parameter dependence of the 
interaction, and in section || we will discuss how the speed of the traveling wave varies with 
impact parameter, averaging out to the velocity observed in the experimental data. 

The paper is organized as follows. In the next two sections we briefly discuss the 
DIS cross section, the concept of geometric scaling, and the dipole cascade approach to 
DIS . In section ^| we discuss the effect of the charm contribution together with the effect of 
saturation on the total cross section. In section |5[ we present the two simple approximations 
to the full model as mentioned above. The scaling properties of the charm structure 
function are studied in section || and in section [7| we concentrate on the region Q 2 < Q 2 
and study the effects of confinement and finite quark masses for small Q 2 below 1 GeV 2 
and the scaling properties of the cross section in this region. In section || we investigate the 
behaviour of the scattering amplitude for different impact parameters and how the scaling 
feature varies with impact parameter. Finally, in section ^, we reach at our conclusions. 



2. DIS and Geometric Scaling 



In the dipole description of DIS the virtual photon, long before the interaction with the 
proton, splits into a qq pair which then interacts with the proton. The transverse separation 
between the quark and the antiquark in such a dipole is denoted r, and their fractions of 
the 7* longitudinal momentum z and 1 — z. The coupling of the 7* to the qq pair is well 
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known and the leading order result reads 
6a„.. 

} 



m 2 r 



\Mz,r)\ 2 = ^P£4QV(1 - z) 2 K 2 (y/z(l-z)Q* + 



+ ml Kl[Jz{l- z)Q 2 + m 2 f r)\ . (2.1) 



Here tpL and denote the longitudinal and transverse wave functions respectively. Kq 
and K\ are modified Bessel functions and the sum runs over all active quarks flavours, 
with mass rrif and electric charge et. The j*p total cross section can then be written as 



o*? p = y d 2 r y ^{|VL(z,r)| 2 + |^T^,r)| 2 }a dp (z,r). (2.2) 
This cross section is related to the F% structure function via the relation 

Q 2 .« 



The factor (Td p (z,r) in ( p^ ) above denotes the dipole-proton cross section. In the GBW 
model this cross section is assumed to have the form given in eq. ( |1 . 1[ ) . If we only consider 
the three light quarks and neglect their masses, we see that the dipole size can be rescaled 
r — ► r/Ro, which implies that the result only depends on the scaling variable r = Q 2 /Q 2 S . 

The charm quark is known to give a significant contribution to the cross section, and 
the heavy charm quark must obviously have a large effect on the scaling properties. This 
will be discussed further in sections § and |(| However, also for the light quarks the finite 
masses will have significant effects for small Q 2 , which will be discussed in section 0. 

3. The dipole cascade model for DIS 

We will in this section shortly describe the model we use to study the effect on geometric 
scaling from different features in the QCD evolution. 

DIS at small x is dominated by gluonic cascades due to the 1/z singularity in the 
gluon splitting function. At large Q 2 the cascade can be described by DGLAP evolution, 
where the gluons are strictly ordered in transverse momentum. For smaller Q 2 also non- 
ordered gluon chains are important, and the fej_-ordered DGLAP evolution is replaced by 
the cc-ordered BFKL evolution. The ordering in k± when Q 2 — > oo and the ordering in 
energy when x — > follow both from an ordering in angle or, equivalently, in rapidity. 
Such an ordering is a consequence of soft gluon interference, and is the basis for the CCFM 
model [18, 19] , which reproduces DGLAP and BFKL evolution in their respective domains 
of applicability with a smooth transition in between. The CCFM model is reformulated 
and generalized in the Linked Dipole Chain (LDC) [20] model. Here some emissions, which 
are treated as initial state radiation in CCFM, are instead included as final state radiation, 
with the result that the "non-Sudakov" formfactors disappear, and the cascade becomes 
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symmetric when exchanging the role of the projectile and the target. The remaining 
gluons (called primary gluons in ref. [20] and backbone gluons in ref. [21]) define totally 
the structure of the final state. 

At very high energies the density of gluons becomes very large, and nonlinear effects 
are needed to tame the exponential growth in the above linear evolutions, which would 
otherwise break the unitarity limit. Saturation and multiple interactions are easier to 
take into account in a formulation in transverse coordinate space since at high energies the 
transverse coordinate does not change between repeated subcollisions, while these collisions 
do change the transverse momenta. This is exploited in the GBW dipole model and the 
dipole cascade model by Mueller [15-17]. 

The evolution in Mueller's model reproduces the leading order (linear) BFKL evolution. 
It starts from a colour singlet qq pair. The quark and the antiquark emit gluons coherently, 
forming a colour dipole. The original dipole is then split in two dipoles formed by the qg and 
gq systems. The new dipoles split repeatedly forming a dipole cascade. In the leading log 
approximation the probability per unit rapidity Y for a dipole with transverse coordinates 
x and y to split emitting a gluon at z is given by (to leading log accuracy Y = In 1/x and 
the true rapidity are equivalent) 



These singularities have to be screened by a cutoff, but the small dipoles have also a small 
probability to interact with the target, and therefore the total cross section is finite when 
the cutoff goes to zero. This implies that a lot of non-interacting virtual dipoles are created 
in the process, which also makes computer simulations difficult [22,23]. 

It is well known that a significant part of next to leading corrections are related to 
energy-momentum conservation [24], and in refs. [25,26] it is demonstrated that energy 
conservation has a large effect on the small x evolution. Relating the dipole size r to 
l/k± there are great similarities between Mueller's model and the LDC model. In ref. 
[13] these similarities were used to implement energy conservation in the dipole cascade 
formalism. The result is a dynamical cutoff for small dipoles, and the remaining emissions 
are ordered in both light cone variables p + and p_, similar to the ordering in the LDC 
model. A dominant fraction of the virtual emissions is here eliminated, leaving mainly 
the "primary gluons" mentioned above. As a result the exponential growth for small x is 
significantly reduced. It also greatly improves the efficiency of the MC simulation, removing 
the difficulties encountered in refs. [22,23]. 

Mueller's cascade includes saturation effects from multiple collisions in the Lorentz 
frame chosen for the calculation, but not the effects of pomeron merging in the cascades. 
The result is therefore not Lorentz frame independent. In ref. [14] we improved our model 
by allowing (colour suppressed) recouplings of the dipole chain during the evolution, a 
"dipole swing", which leads to an almost frame independent formalism. In this paper we 
also introduced a simple model for the proton in terms of three dipoles. 





The splitting probability in eq. Q3.1| ) is singular for small dipole sizes x — z or z — y. 
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Figure 1: Our full results for the j*p total cross section plotted as a function of the scaling 
variable r = Q 2 /Q s GW . Here saturation effects are included both in the collision between the 7* 
and the proton cascades and within the evolution in each individual cascade, via the so called dipole 
swing. We have used a charm mass of 1.4 GeV and an effective light quark mass of 60 MeV for the 
u, d and s quarks. Data points are taken from [27] and [28]. Results are presented for Q 2 ranging 
from 0.3 GeV 2 to 90 GeV 2 , and with the same kinematics for the data and the model. 

The dipole splitting is calculated in perturbative QCD, and therefore the model is 
meant to work in the perturbative regime, which means not too small Q 2 . In figure [j] 
we show results for a 7 p as a function of the scaling variable r = Q 2 /Q 2 for different Q 2 
above 0.3 GeV 2 . The results include both the dipole swing and multiple collisions and for 
Q s {x) we here use the definition by Golec-Biernat and Wusthoff, given by (|1.2| ). For these 
Q 2 - values a scale-breaking effect is obtained from the charm mass, for which we use the 
value 1.4 GeV. The theoretical result is presented for the same kinematical variables as the 
experimental data, and we see that there is a very good agreement between theory and 
data. For small Q 2 the result is sensitive to confinement effects and effective quark masses. 
These problems will be further discussed in sec. [/]. 

4. Effects of saturation and the charm contribution 

The effects of saturation and of the large charm quark mass are illustrated in figure The 
results correspond to Q 2 between 0.75 and 90 GeV 2 , and therefore the light quark masses 
can be neglected. The dotted lines show the results of the full model, presented in fig. |l|. 
The solid lines show the result when the charm mass is set to zero. The deviation between 
these curves therefore shows the effect of the charm quark mass. Neglecting also the swing 
gives the long-dashed curves, and finally including only the single pomeron term in the 
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Figure 2: The top three curves show the j*p total cross section obtained from our Monte Carlo 
simulation for Q 2 ranging from 0.75GeV 2 to 90GeV 2 and for zero quark mass for 4 flavours. The 
solid lines show the results for evolution including the dipole swing and also multiple collisions, the 
long-dashed lines show the results without the swing but including multiple collisions, while the 
short-dashed lines are results without the swing and also without multiple collisions. It is seen that 
all three results satisfy geometric scaling. The bottom curve represents our full results shown in fig 

B 



collisions gives the short-dashed curves. The difference between these and the solid lines 
therefore represents the effect of saturation. 

We first note that saturation effects from multiple collisions and the swing have a 
relatively small effect for r > 1, but grow for smaller r, and reduce the cross section by 
approximately a factor 2.5 for r = 0.1. The effect of the charm quark mass is, as expected, 
also largest for smaller r. The charm contribution is about 30% for r = 100 and 10% for 
T = 1, which should be compared with 40% for zero mass charm quarks. 

We note in particular that also the one-pomeron result satisfies geometric scaling, even 
down to t = 0.1. Thus the scaling feature alone is not enough to prove the existence of 
saturation. The scale-breaking effect from the charm contribution will be further discussed 
in sec. ^. The dominance of the light quark contributions, and the strong correlation 
between Q 2 and x due to the limited HERA energy, imply that this effect is not so evident 
in the results shown in figs. |l| and although we can see a small shift between the curves 
for different Q 2 in fig. ||. 

Although saturation reduces the result for r < 0.5, there is a more clear break in the 
experimental data. We note, however, that the experimental points, which show this break, 
all correspond to small Q 2 -values below 0.5 GeV 2 . Here we expect also non-perturbative 
effects to be important. The wave functions in eq. ( [2.1| ) extend to very large transverse 
separations r, and such large dipoles must be suppressed by confinement effects. We also 



-8- 



note that energy limitations imply that for fixed Q 2 the experimental points lie in a rather 
small x-interval, which implies that the points for different Q 2 have a limited overlap. The 
scaling features for small Q 2 will be further discussed in sec. [7]. 



5. Understanding geometric scaling in the linear evolution 

The saturation model gives a motivation for geometric scaling in the kinematic range 
Q 2 < Q 2 , but it does not give a reason for the observed scaling behaviour for Q 2 > Q 2 
at experimentally feasible energies at current accelerators. The arguments in ref. [6] imply 
that the cross section should scale as a function of Q 2 /Q 2 (x) within a wide kinematic 
region, 

i < HQ 2 /Ql) « HQ 2 /A 2 QCD ), (5.1) 

which is dominated by BFKL diffusion. Here Q 2 (x) is defined as the scale where the 
scattering probability is of order 1. As mentioned in the introduction this is a consequence 
of the fast growth towards small x and wide diffusion to large k± in the leading log BFKL 
evolution. Both these effects are strongly damped by non-leading effects [7], which therefore 
reduce the region dominated by linear BFKL evolution. (We also note that both references 
[7] and [6] point out that with a running coupling the BFKL evolution leads to a saturation 
scale which satisfies InQ 2 ~ \/Y rather than lnQ^ ~ Y.) For larger Q 2 , fcj_-ordered 
evolution should be dominant and Kwiecihski and Stasto [8] have assumed that DGLAP 
evolution is applicable in the whole region Q 2 > Q 2 , with scaling initial conditions at Q 2 = 
Q 2 . They then found that although not perfectly scaling, the result showed approximate 
scaling in the same kinematical domain specified by eq. ( |5.1[ ). Their arguments do, however, 
not explain the shape of the saturation line. 

We will below argue that scaling is expected both in the region dominated by BFKL 
diffusion and in the /cj_-ordered double leading log (DLL) regime at larger Q 2 . (In the 
DLL regime this is not obvious from the analytic expressions, but follows from a numerical 
analysis.) Furthermore we find that the transition between the DLL (/cj_-ordered) and 
BFKL (/c_i_-non-ordered) regimes is given by Q 2 = Q 2 imit oc x Xbfkl , where Xbfkl is the 
exponent in the solution to the BFKL evolution, estimated to be around 0.3. This implies 
that Qf imit is very close to Q 2 GW , leaving only a very small region where the linear BFKL 
evolution is dominating. 

The qualitative features of the QCD evolution at small x are present already in the 
leading log 1/x results, and we will in sec. |5.1| discuss a toy model describing the BFKL- 
DLL transition [29]. Non- leading corrections from e.g. energy conservation and a running 



coupling are important for the quantitative result. In sec. |5.2| we try to isolate the most 
important features of the full dipole cascade MC, to get an intuitive picture of geometric 
scaling in the non-saturated region. 

5.1 Leading log approximation 

When both 1/x and Q 2 are large the gluon distribution is given by ladders which are ordered 
in k± and where the splitting function is dominated by the 1/z pole. This corresponds to 
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Figure 3: The cross section er 7 » p ~ g/Q 2 with the gluon density g given by equation for 
different Q 2 , varying from Q 2 = lOOGcV 2 to Q 2 = 1.5 • 10 10 GeV 2 , and plotted as a function of 
Q 2 /Q 2 , with Q s parametrized according to fll.2|), but with A equal to 0.7. 



the double log approximation, and for a running coupling the gluon density is given by 



g(x,Q 2 ) ~ exp (2\/ao lnlnQ 2 In T/x ) , where a(Q 2 ) = 



ln(Q 2 /A 2 )" 



(5.2) 



This expression does not scale exactly with Q 2 oc x~\ Neglecting the very slow variation 
of In In Q 2 the cross section ~ g/Q 2 scales with Q 2 oc exp(A'-v/ln 1/x), with some parameter 
A', and in ref. [30] it is pointed out that the experimental data can be equally well fitted 
with both these expressions for Q s . This is also seen in fig. ||, which shows the cross section 
ijT'p ~ g/Q 2 vs. t with g from eq. (5.2) and r defined from eq. ( |1.2j ) but with A = 0.7. 

For limited Q 2 we are instead in the BFKL regime dominated by /c^-non-ordered 
chains, where the gluon density is growing as a power l/x x for small x, with A of the order 
0.3. Multiplying by 1/Q 2 therefore gives directly the scaling cross section 

afP ~ ^| = r- 1 (5.3) 

with Q 2 S oc l/x - 3 . 

We now also want to argue that, as described in ref. [29], the line Q 2 = Q 2 imit corre- 
sponding to the separation between ordered (DGLAP-like) chains and unordered (BFKL- 
like) chains, is also close to Q 2 GW - In the double log approximation the (non-integrated) 
gluon distribution G is given by the fc^-ordered expression 



4a, dk 2 



-L,i 



3vr k 2 



6(k A 



-i)—0(x i 



x/) I 5(x - x n )5(k]_ - k 2 L n ). 

(5.4) 
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With the notation k = ln(l/xj) and Ki = \n(k\Jk 2 ) (and k = ln(fc^/A 2 ), I = ln(l/x)) we 
get for a fixed coupling the result 

( n ~ K n Anl/x ~\ 

G ~E|ny adKiO(Ki - Ki-i) ■ Yl J dU0(li-U-i)> 

E K n in , . 

a n - — ■— = / (2 v / alnQ 2 lnl/x) ~exp(2 v / alng 2 lnl/x . (5.5) 
n\ n\ V / 

n 

For a running a s we have instead of dk\ Jk\ i = dni a factor dni/ Ki = dln/^, which then 
gives the result in eq. ( |5.2j ). 

In the BFKL region with small x but not so large Q 2 , the fc^-ordered phase space in 
eq. ( [5.4D becomes small, and chains which are not ordered in k± give important contribu- 
tions. The BFKL evolution can be formulated in different ways. Expressed in terms of 
the primary [20] or backbone [21] gluons a step downwards in k± is suppressed by a factor 
k±i/k±i-i [20,31]. We note that this implies that a maximum A:^-value in the chain will 
contain the factor dk 2 ^^/ k\ which can be interpreted as a hard parton-parton subcol- 
lision with the expected cross section proportional to dt/t 2 . Expressed in the logarithmic 
variable k, a step down is consequently suppressed by a factor exp(Kj — = exp(— 8k). 
This implies that the effective range allowed for downward steps corresponds to approxi- 
mately one unit in k. Consequently we find that the phase space limits K{ ^ in eq. (|5.5| ) 
is replaced by [29,31] 

Ki ^ Ki-i — 1. (5-6) 



For a fixed a s the transverse momentum integrals giving K n /n\ in eq. ( |5.5|) will be 
replaced by 

f fr d^d^ - k^ - 1) « impl (5.7) 



When k is very large we recover the DLL result in eq. ( p.5|) , but for smaller values of k we 
find instead using Stirling's formula 



which implies 



with 



(k + n) n n n , 
k small => ~ e n 5.8) 

TV. TV. 



G [aeln(l/x)] n _ ^ainfi/s) _ 1 ,g q n 



A = ea« 2.72a. (5.10) 

In this range the chain corresponds to a random walk in Ink 2 ^. The result should be 
compared with the result from the leading order BFKL equation, which gives 

A = 41n2a w 2.77a. (5.11) 

We see that this simple picture describes the essential features of BFKL evolution. 
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We note here in particular that the boundary between the domains dominated by 
fcj_-ordered chains (eq. (|5.5|)) and k± -non-ordered chains (eq. fl5,8j )) is determined by the 
relation 

ln *aiimit = edln(l/x). (5.12) 

We see that if we replace the LL A-value by a lower value ~ 0.3, as indicated by the 
experimental data, then fc5_ limit is similar to Q 2 GW (x) from the early fit by Golec-Biernat 
and Wusthoff [2,3]. This shows that the line Q 2 = Q 2 qw(x) (corresponding to r = 1) is 
actually close to Qf imit , which represents the separation between /cj_-ordered (DGLAP-like) 
chains and non-ordered (BFKL-like) chains. 

5.2 The dipole cascade 

The qualitative features in the LLA in sec. |5.1| are modified by energy conservation and 
other non-leading effects. We will here study how the scaling behaviour can be seen in the 
dipole cascade model presented in sec. ^, trying now to isolate the most important features 
of the full MC simulation. 

Let us look at a dipole with size r, which after a rapidity interval Ay splits in two 
dipoles with sizes r> and r < , with r > larger than r< . Fig. [| shows the MC results for the 
average values of the ratios r>/r and r < /r for different values for the "time", Ay, it takes 
from the formation of the dipole till it splits. The probability distribution dP/dAy is also 
shown. We see here that a dipole typically splits in two dipoles, where one has the same 
size as the parent while the other has just half its size. This result is independent of the 
step Ay and of the photon virtuality Q 2 . This also shows that it is independent of the size 
of the parent dipole. We note that these results depend crucially on the energy-momentum 
conservation in the evolution. Without this constraint the typical splitting would occur 
within a very small rapidity step giving one virtual dipole with a size close to the necessary 
cutoff. When the cutoff goes to zero both Ay and r < /r would also go to zero. Essential 
for the simple result in fig. |I| is the conservation of both p + and p_. This is seen in fig. || 
which shows the corresponding result obtained with only p + -conservation. 

Let us then assume that each dipole r splits after a typical Ay, which from the distri- 
bution dP/dAy is found to be around 1.8 units, into two dipoles with sizes r and r/a, with 
the parameter a of the order of 2. If Y denotes the total rapidity range the number of steps 
will be A = Y I Ay and the number of dipoles 2^ = 2^ Y /^ y \ Starting from an initial dipole 
ro the number of dipoles having size tq/oJ 1 will be equal to ( ) , with n = 0, 1, . . . , N. 

To study the scaling behaviour of the 7*p cross section within this approximation, 
we can treat the photon as a dipole with size r~ = 1/Q. In the linear region the initial 
proton can also be treated as a single dipole with size r p = 1/A. The cross section for the 
scattering of two dipoles, r\ and r 2 , is given by [17] 

CTdd(ri,r 2 ) = 2ira 2 s r 2 min {l + ln(^^)} (5.13) 

where r max (r m i n ) is the largest (smallest) of the two colliding dipoles r\ and r 2 . We study 
the collision in a frame where only the proton is evolved, and the photon therefore treated 
as a single dipole. 
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_Q 2 =1 GeV 2 
Q =100 GeV 2 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 

Ay 

Figure 4: The top figures show how the average values of the ratios r</ro (left) and r>/ro 
(right) for splitting events vary with the rapidity separation Ay. It is seen that both these ratios 
are approximately constant independent of Ay, and that r>/ro is ~ 1. The bottom figure presents 
the Ay-distribution dP/dAy. This distribution has a peak around Ay = 1 and falls off exponentially 
for larger Ay-values. 




0.1 ' 1 1 1 1 1 1 0.1 

0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 

Ay Ay 

Figure 5: The averaged ratios (r</ro) and (r>/ro) as functions of Ay in a simulation where 
P- is not conserved in the evolution. We see that ^-conservation is not enough to get the almost 
constant ratios seen in figure ||. 

When Q is larger than a N -A the photon dipole is smaller than all dipoles in the proton 
cascade. The resulting cross section is then given by 

,>v>_^fm n -.l Q 



n=0 



™;Q- 2 2 "l^- Nina + ln^ } oc |f jl + \ ln(^-) \. (5.11) 
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where in the last expression we have used a = 2, iV = Y/Ay and Q 2 C = A 2 x~ ln2 ^ Ay ~ 
A 2 a; -0 ' 38 . The result is obviously scaling, and the exponent 0.38 is not far from the ex- 
perimental fit around 0.3. We note also that in this case the dominating contribution 
comes from dipole chains where the dipoles are ordered in size and where the last dipole 
in the cascade is larger than the photon dipole. This just corresponds to the dominance of 
fc_i_-ordered DGLAP-type ladders. 

For smaller Q- values the curly bracket in eq. ( 5.14| ) overestimate the contribution from 



small proton dipoles represented by n-values for which a n -A > Q. The largest terms in the 
sum are obtained when the binomial factor has its maximum, i.e. for n ~ N/2. Therefore 
the result in eq. ( |5.14 ) is a good approximation as long as Q > a N ' 2 -A. Indeed, using 



Stirling's formula we can write (^y 2 ) ~ ^ N ■> an d in the saddlepoint approximation we again 
arrive at the result given in eq. ( |5.14D . We note that the constraint Q > a N / 2 -A corresponds 
just to Q 2 > Q 2 C , which thus is the limit for the dominance of fe^-ordered, DGLAP-type 
evolution chains. 

Below this region, i.e. for Q 2 < Q 2 C , it is not a good approximation to neglect the 
contributions from scatterings where the proton dipoles are smaller than the photon dipole. 
The full expression, including all contributions, can be written as 



n=0 v J K J n=m+l v 7 ^ J 



(5.15) 



where m = — • This expression is more complicated, but we show in figure || the result 
of a numerical evaluation expressed with a scale Q 2 C oc x~ 0A . We see that scaling is indeed 
satisfied for a large range of values for Q 2 , also in the kinematic region dominated by chains 
which are not well ordered in dipole size or in transverse momentum. 

We note that our toy model has important similarities with the empirical model pre- 
sented in ref. [32]. In this paper it is demonstrated that geometric scaling follows from an 
assumption that the dipole cascade in the proton is dominated by dipoles with size Qj 1 
(while the virtual photon can be represented by a dipole with size Q^ 1 ). Our toy model 
approximation of the full MC cascade has just this feature; the dipole multiplicity is given 
by the binomial coefficient (^), with a maximum for n = N/2 corresponding to r = 1/Q SC . 



The result in eq. ( 5.14 ), which is obtained for Q > Q sc using the saddlepoint approxima- 
tion, is thus identical to the corresponding result in [32]. Ref. [32] also points out that 
if events with Q < Q s can be described by an analogous cascade evolution of the initial 
photon dipole, then the expression Q/Q s a' y * p (Q/Q s (x)) is symmetric under the exchange 
Q <-> Qs(x). This symmetry was initially observed in the HERA data by ref. [1]. For these 
smaller Q- values (Q < Q sc ) the dominant contributions from n ~ A^/2 are contained in the 
second term in eq. ( |5.15| ). Using again the saddlepoint approximation we also obtain the 
symmetric result 

a^P(Q 2 ,Y) « 27ra 2 2 N a- N A- 2 jl + m(^^) j « 2™^ ji + m(^)j . (5.16) 
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Figure 6: The toy model cross section in flSUSl ) plotted as a function of the scaling variable with 
Q 2 C = A 2 -x~ aA , with A 2 = O.lGeV 2 , and for different Q 2 . The normalization here is not of interest 
and all prefactors has been simply put to 1. We indeed see that the result shows scaling, for low 
and high Q 2 alike. 



In conclusion we find that Q 2 = Q 2 SC ocx , with A approximately equal to the BFKL 
exponent, specifies the limit between dominance of /cj_-ordered and k± -non-ordered chains, 
and that a simple toy model having this property gives the qualitative features of the scaling 
dynamics for both large and small r. We also see that the toy model approximation to 
the full MC simulation of the dipole cascade model have important similarities with the 
empirical model in ref. [32], and gives a scaling result with a scaling exponent not far from 
what is observed at HERA. 



6. Scaling features in the charm contribution 

It is well known from HERA data that charm quarks contribute a significant part to the 
total cross section. As discussed in sec. ||, the large charm quark mass modifies the scaling 
properties, and it is seen that the HERA charm data do indeed not scale as a function 
of r = Q 2 /Q 2 GW (x). As pointed out in ref. [33] they do, however, scale quite well as a 
function of the modified scaling variable r c = (Q 2 + Am 2 )/Q 2 

From the photon wave function in ( |2.1| ) we see that the term proportional to the Bessel 
function K 2 in \4>t\ 2 only contains Q 2 and mj in the combination z(l — z)Q 2 + m 2 . If 
z- values around 1/2 dominate, we would expect that the charm contribution is gradually 
switched off when Q 2 is of the order Am 2 or smaller, which necessarily leads to a breaking 
of geometric scaling. The sum of the terms proportional to Kq in \4>t\ 2 and \iPl\ 2 is 
proportional to [42 (1 — -z)]z(l — z)Q 2 +m 2 f , and for z close to 1/2 the square bracket equals 
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Figure 7: The total charm cross section, Oj, + a\ plotted as a function of r c defined in (6.1) with 
A = 0.35 and m c = 1.4 GcV. Results are shown for Q 2 between 2 and 130 GeV 2 . We see that the 
result scales fairly well with this scaling parameter. Data points are taken from ref [34]. 



1, and we get the same factor z(l — z)Q 2 + nil as before. From these features we may 
expect that the charm contribution scales approximately with a modified scaling parameter 
where, for example, we replace Q 2 by Q 2 + n-m 2 in the definition of r, for some number n 
which should be close to or a little larger than 4. 

In fig. [7| we show data for the charm cross section from Zeus [34] and our MC obtained 
for m c = 1.4 GeV. In this figure we have used the scaling variable 

t c = (Q 2 + 6m 2 c )/Ql GW {x) with x = (Q 2 + 6m 2 c )/W 2 , (6.1) 

obtained with n = 6 in the definition of r c , which we find gives somewhat better scaling 
properties. This implies that r c = (1 + Qm 2 / 'Q 2 ) 1+x -tqw ■ We here first note that the MC 
agrees quite well with the data, being only a little high for the largest r-values. Secondly 
we confirm the observation in ref. [33] that the charm cross section does scale well as a 
function of such a modified scaling variable. 

The value n = 6 chosen for the scaling variable in eq. (|6.1| ) is not crucial for the scaling 
behaviour, which works rather well also for n = 4. We note, however, that replacing xby x 
in Q s only contributes a minor part to the difference between r c and tqw-, represented by 
the factor (Q 2 + 6m 2 ) x . More important is the replacement Q 2 — > Q 2 + 6m 2 in eq. (|6.1|). 
This is illustrated in fig. [8|, which shows that scaling does not hold if we replace r by 

We want here to point out that the effect of the charm mass is reduced very slowly for 



-j Q ' 1 1 1 1 ' ' ' ' I 1 1 1 1 ' ' ' ' I 1 1 1 1 ' ' ' ' 
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x GW (xbar) 

Figure 8: Here we plot the charm cross section as a function of Q 2 /Q s (x) with the same param- 
eters as in fig. 0. In this case we see that the result does not scale, which shows the importance of 
replacing Q 2 with Q 2 + n ■ m 2 , n ~ 4, as argued in the text. 

larger Q 2 . Integrating over the dipole size, r, we obtain from a dimensional analysis 

J d 2 rK 2 (er)a(z,r) oc i , e = yjz(l - z)Q 2 + mj . (6.2) 

When 4Q 2 is small compared to m 2 we have e ~ mj, and the charm cross section ought 
to scale as 

cr T oc — y, a L oc —j. (6.3) 

m L c mil 

We see here that the longitudinal contribution depends on two separate scales Q 2 and m 2 . 
In the other limit, when z(l — z)Q 2 > m 2 we can neglect m c . However, since z(l — z) 
can take on arbitrarily small values for fixed Q 2 , we the effect of the charm mass can be 
important also for high Q 2 . In the MC results presented in fig. [7| the reduction due to 
the charm mass is about 30% for Q 2 = 90GeV 2 and 55% for Q 2 = 15GeV 2 . This slow 
decrease of the mass effect is also seen in fig. |2[ (We should here also remark that the 
transverse wave function tpx is not normalizable due to the singularity K 2 (r) ~ 1/r 2 at 
small r. However, at small r the dipole-proton cross section behaves like r 2 and thus the 
result is still finite.) 

7. Interaction at smaller Q 2 

Due to the limited energy in the HERA accelerator the experimental data in the region 
r < 0.5 are all obtained for rather small virtualities Q 2 , where non-perturbative effects 
must be expected. 
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Figure 9: The j*p total cross section plotted as a function of the scaling variable r for low Q 2 
values. Points are data from the HI [27] and ZEUS [28] collaborations, while the lines are results 
obtained from our Monte Carlo. We can here see a successive suppression of the cross section for 
smaller Q 2 as a result of the finite quark masses. 



7.1 Can the perturbative dipole formalism be used for Q 2 below 1 GeV 2 ? 

For small Q 2 the wave functions in eq. (2.1) extend to very large transverse separations r, 
and such large dipoles must be suppressed by confinement effects. The interaction of these 
photons is usually described in terms of two separate components, a vector dominance 
contribution and a direct coupling to the quarks. We will here test if it is possible to 
represent the interaction of these photons with a finite mass for the light quarks, which 
effectively suppresses the contribution from the very large dipoles. Because the contribution 
from the strange quark is suppressed by its smaller electric charge, and therefore relatively 
small, it is not possible to study in any detail the effect of the strange quark mass. We 
therefore use a single quark mass, mi, for all the light quarks, u, d, and s. The results in 
fig. H for Q 2 in the range 0.3-3.5 GeV 2 are obtained for mi = 60MeV and m c = 1.4GeV, 
and we here show the model results for the same combinations of x- and Q 2 -values as in 
the experimental data. We see that the experimental data are very well reproduced by 
the model calculations, which gives some support to the application of this perturbative 
description also for these small virtualities. 

The effective quark mass is quite small, and we would have expected a larger suppres- 
sion for small Q 2 . The value 60 MeV is smaller than Aqqd, and also smaller than the 
effective mass obtained in ref. [35] in an analysis of the vector-current two-point function. 
Ref. [35] studies a model with an effective quark mass which varies with Q 2 , becoming 
smaller when Q 2 is increased. The agreement in fig. ^ would actually be improved by such 
a varying effective mass, but we do not believe that the accuracy of our model is sufficient 
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Figure 10: In this figure we extend the result for Q 2 = 2GeV 2 to smaller r-values. We can 
here see a deviation from the scaling behaviour at higher energies, about a factor 1.4 for r « 0.07 
corresponding to x ~ 3 • 10~ 9 for Q 2 = 2GeV 2 . In this plot the deviation from the scaling curve is 
to a large extent due to the effective light quark mass. 



to claim that this improvement is significant. We have therefore here only used a constant 
effective quark mass in the wave functions in (2.1). 



7.2 Is geometric scaling obeyed for Q 2 /Q 2 S < 1? 

The effect of a finite quark mass is approximately a multiplicative factor which suppresses 
the cross section for smaller Q 2 . As discussed in sec. |6| the result of the quark mass in the 
photon wave functions corresponds roughly to a suppression by a factor Q 2 /(Q 2 + 4m 2 ). 
The fact that the model results in fig. |^ look as a scaling function is therefore a consequence 
of the strong correlation between x and Q 2 in the experimental data. The finite energy 
in the HERA accelerator implies that for small values of r there are only data for (x, 
Q 2 )-values within a very small interval. With a future accelerator with higher energy, one 
could also reach smaller T-values keeping Q 2 > lGeV 2 . For these Q 2 we do not get much 



suppression from the light quark masses, and fig. |10| shows that the scaling curve will lie 
somewhat above the present HERA results. In the figure we show the extend the curve for 
Q 2 =2GeV 2 to t ~ 0.07. As can be seen, the difference is increasing for smaller r-values 
and is about a factor 1.4 for r rj 0.07. We also expect to see the scalebreaking effects of 
the charm mass for lGeV 2 < Q 2 < lOGeV 2 , while for higher Q 2 these effects should be 
gradually switched off. This cannot, however, be observed at accelerators in the foreseeable 
future. 

As a conclusion of this section we predict that at higher energies the data will scale 
at larger values for the cross section than the present HERA data for r < 1. We should, 
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however, point out that the analysis in this section only holds for the particular value 
A = 0.3 for the saturation scale in eq. ( |1,2| ). For a different A the scaling behaviour will be 
different and it turns out that the deviation from the scaling behaviour would be somewhat 
reduced with a higher A. 



8. Impact parameter dependence of the traveling wave 

In this section we will look at the impact parameter dependence of the dipole-proton 
scattering amplitude T dp (b), defined as T = 1 — S. We expect that the contribution from 
small dipoles is reduced for large impact parameters, which may result in different scaling 
behaviour for central and peripheral collisions. This feature may be of special interest in 
studies of diffraction. 

The amplitude T dp is related to the total dipole-proton cross section by the relation 



a dp (z,r) 



2 J d 2 bT dp (b). (8.1) 



To get the amplitude for photon-proton scattering, the amplitude T dp (b) has to be weighted 
by the photon wave function in eq. (|2.l|). The average value (T) will then depend on the 
virtuality Q 2 and the impact parameter b. 

In the introduction we mentioned the similarity between the high energy evolution 
equations for T and a certain type of equation (or rather a class of equations) from statistical 
physics, known as the FKPP equation. Neglecting the impact parameter dependence, the 
BK equation for the amplitude T{l,Y) (with I = Ink 2 ) is analogous to the FKPP equation 
for the function u(x,t), with the identifications I <-> x and aY «-> t. The asymptotic 
solution, as t — > oo, u as (x,t) = u(x — vt), then corresponds to T as {l,Y) = T(l — XY) 
in QCD. If we define Q 2 = exp(Al") as before the solution satisfies the geometric scaling 
relation T(l, Y) as = T(lnk 2 /Q 2 ). There are a number of conditions which must be satisfied 
in order to obtain an asymptotic solution of this form, and for a short review on traveling 
wave solutions in QCD we refer to [36]. 

As our model satisfies geometric scaling we can expect that the function (T)(lnQ 2 ) 
plotted for different Y- values looks like a traveling wave. This can be seen in fig. |ll| which 
shows results for Y = 6,8, 10 and 12. In this figure we also show the result for different 
impact parameters b. If we look at different points with the same (T)-value, the velocity 
is determined by the relation 

fh\ AlnQ2 

v(b) 



AY 



The results in fig. |ll] correspond to the following velocities when saturation effects are 
included in the amplitude: v(b = 0) = 0.28, v(b = 0.5/m) = 0.35 and v(b = 1/m) = 0.37. 
(For the one pomeron amplitude we instead obtain v(b = 0) = 0.41, v(b = 0.5/m) = 0.46 
and v(b = 1/m) = 0.48) Thus the velocity varies significantly with the impact parameter, 
being smaller for central collisions and larger for peripheral collisions. As the scattering 
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Figure 11: The average scattering amplitude as a function of ln(Q 2 /iGeV 2 ) f or Y = 6, 8, 10 and 
12 and impact parameters 6 = 0, 0.5, lfm. Both the amplitude containing full saturation effects and 
the linear one pomeron amplitude are shown. For the amplitude including saturation we obtain 
from here the velocities v(b = 0) = 0.28, v(b = 0.5/m) = 0.35 and v(b = lfm) = 0.37. 



probability is largest for the central collisions, these results seem to be reasonably consis- 
tent with weighted average over impact parameters corresponding to the velocity v ~ 0.3 
observed for the total cross section. 

In fig. |ll| we can see that also the one pomeron amplitude, which contains no saturation 
effects, seems to exhibit the form of a traveling wave. We also see that we need quite high 
energies, corresponding to x ~ 6 • 10~ 6 for Q 2 ~ 4GeV 2 and 6 = 0, before the one pomeron 
amplitude reaches the unitarity limit T = 1. However, it is also seen that saturation effects 
reduce the amplitude by about 10% for larger Q 2 and about 40% for small Q 2 ~ lGeV 2 
already at x ~ 10~ 3 for 6 = 0. This result is consistent with the previous results presented 
in fig. |. 



9. Conclusions 

DIS data from HERA show a striking regularity as <7 7 * p is a function of the ratio r = 
Q 2 /Q 2 (x) only, with Q 2 (x) given by eq. ( |E^ ) [1-3]. Such a geometric scaling has been 
expected in the range r < 1, as a natural consequence of saturation when the gluon 
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density becomes very high, and there has been a lot of discussion in the literature whether 
saturation has or has not been observed at HERA. Modifications of the saturation model 
including DGLAP evolution [4] improves the agreement with data for larger Q 2 , but it 
has not been equally obvious how geometric scaling follows in a natural way from QCD 
evolution in the non-saturated domain. The traveling wave solutions to the nonlinear 
evolution equations in QCD do predict geometric scaling also for r > 1 [9,10], but these 
solutions are valid only at extremely high energies, far beyond what is available at the 
HERA accelerator. 

In this paper we have tried to shed some light on these questions: What is the reason 
for scaling for r > 1, and is saturation indeed the reason for geometric scaling for r < 1? 
Supporting the saturation idea is the fact that the scaling curve seems to have a break 
just around r = 1. There are, however, also other effects which contribute to a break in 
this region. In ref. [29] it is shown that the transition between ^-ordered (DGLAP-like) 
chains and fc^-non-ordered (BFKL-like) chains is expected close to Q 2 . Secondly the finite 
energy in the HERA machine implies that experimental data for r < 1 are only available 
for Q 2 < lGeV 2 , where nonperturbative effects begin to be important. 

To study the influence on the scaling from different dynamical effects we have used 
the dipole cascade model presented in ref. [14]. This model is based on Mueller's dipole 
cascade model, but includes energy-momentum conservation and saturation effects, not 
only from multiple sub collisions but also from pomeron mergings in the cascade evolution, 
and it does successfully reproduce the data from HERA. 

Our conclusion is that scaling is a natural consequence of the dipole evolution in the 
linear region r > 1. Indeed, neglecting saturation the linear evolution exhibits geometric 
scaling at high energies for r- values both larger and smaller than 1. The change from the 
DGLAP to the BFKL regime causes a change in the slope of the scaling curve, but this 
change is rather smooth, without a sharp break at r = 1. The break seen in the HERA data 
is to a significant part caused by saturation, but to an even larger extent by scalebreaking 
effects at low Q 2 . The latter are partly due to the large c-quark mass and partly due to 
non-perturbative effects related to confinement for the light u-, d-, and s-quarks. Therefore 
we predict that the results from a future higher energy machine will show deviations from 
the scaling behaviour in the variable tqw- These results are expected to lie above the 
HERA results for tqw < 1> an d be represented by a curve which is more smooth around 
t~GW = 1- (Recently it has also been predicted that a new type of scaling, called diffusive 
scaling, will occur at very high energies, but we have in this paper concentrated on energies 
which might be within reach at some future accelerator.) 

We study the charm contribution separately, and compare with data for the charm 
cross section. The charm contribution does not scale as the total cross section, but we 
see that it does scale rather well as a function of the variable (Q 2 + 6m 2 )/Q 2 (x). Finally 
we have studied how the scaling feature varies with impact parameter. We here see that 
the effective power A in (|1.2| ) is somewhat smaller for central collisions and larger for 
peripheral collisions, averaging out to the observed value around 0.3. Such a variation may 
be of interest in studies of diffractive scattering. 
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